Hysteresis loops of Co-Pt perpendicular magnetic multilayers 
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We develop a phenomenological model to study magnetic hysteresis in two samples designed as 
possible perpendicular recording media. A stochastic cellular automata model captures cooperative 
behavior in the nucleation of magnetic domains. We show how this simple model turns broad 
hysteresis loops into loops with sharp drops like those observed in these samples, and explains their 
unusual features. We also present, and experimentally verify, predictions of this model, and suggest 
how insights from this model may apply more generally. 
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Magnetic thin films with magnetization oriented per- 
pendicular to the plane of the film have become attrac- 
tive as magnetic bits get smaller and smaller. Cobalt- 
platinum multilayers are thin films with several desirable 
characteristics for future recording media, such as large 
ratio of coercive field (field where M = 0) to saturation 
magnetization and sharp initial drops in their hystere- 
sis loops. Although grain sizes are too large at present 
for applications, they are nevertheless interesting model 
systems, both as high coercivity thin films Q, and as a 
testing ground for Ising model based simulations of hys- 
teresis in real materials. 

These films are also of fundamental interest as interact- 
ing dynamical systems. In particular, the precise onset of 
magnetization reversal, which has been observed before 
in Co-Pt multilayers 0, 0, El 0, Q , suggests the shape 
of the hysteresis loop may be dominated by nucleation. 
This is the case for some magnetic thin films of interest 
in recording [3- C3- IZ3| - Nucleation phenomena are of broad 
interest, e.g. in phase transitions 0, materials synthe- 
sis || , and protein aggregation Q . Unresolved issues for 
these Co-Pt multilayers include the origin of the broad 
tail in the hysteresis loop near magnetic saturation and 
the drop in first-order reversal curves following reversal 
of the magnetic field (Fig. 0. Additionally, differences 
between nominally similar samples are not well under- 
stood. 

Here we obtain the major hysteresis loop and first- 
order reversal curves (measured by decreasing the mag- 
netic field along the major hysteresis loop until a pre- 
defined reversal field, then tracing out a reversal curve 
as the field is again increased to saturation) for two Co- 
Pt multilayer samples, then develop a lattice-based the- 
oretical model to explain key features of these measure- 
ments. The model makes testable predictions concerning 
time-dependence of the reversal curves, which we test, 
and temperature-dependence of the hysteresis loop. Our 
model is remarkable not so much in terms of quantitative 
agreement with experiment, but in the physical insight 
it provides into a variety of experimental results. 

Two Co-Pt multilayer samples were provided by B. 
Terris at IBM. These samples are similar in composition 
to those described by Weller et al. and consist of a 
20 nm Pt buffer layer, followed by 10 Co(0.6 nm)/Pt(l 
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FIG. 1: Major hysteresis loops (envelope) and first order re- 
versal curves for two samples provided by B. Terris at IBM. 
Sample (a) shows sudden sharp drop in major hysteresis loop, 
while sample (b) shows a similar drop as the magnetization 
begins to fall, but then has a broad tail. 



nm) bilayers. Growth temperatures for the multilayers 
were 157°C and 320°C for the samples of Fig. QJa) and 
(b), respectively. Reversal curves (Fig. were taken 
with an Alternating Gradient Magnetometer (Princeton 
Measurements Corporation). The general procedure for 
each reversal curve was to decrease the magnetic field 
every 0.1s until reaching a predefined reversal field. Data 
was then collected at steps along the reversal curve with 
an averaging time at each data point of 0.13s. Variations 
on this procedure were used to test predictions of our 
model and will be discussed below. 

The experimental hysteresis loops and reversal curves 
for these samples exhibit several unusual features. To 
our knowledge, a simple explanation of these has not yet 
been provided. Particularly, we note the broad tail (as 
the magnetization approaches saturation) shown in Fig. 
Gib). Modeling in terms of exchange coupling can pro- 
duce sharp drops due to an energy barrier associated with 
nucleation of reverse domains, but the rounding at the 
tail is typically comparable to the rounding at the onset 
(as the magnetization begins to fall) jl0llll|. We notice 
here that the tail is much broader than the onset. This 
has been observed previously I3-O10] and described as an 
asymmetry around the coercive field . Additionally, we 
notice a drop in the reversal curves in Fig. QJa) following 
reversal of the magnetic field. That is, the magnetization 
continues to fall even after the magnetic field begins to 
increase. Reversal curves are less commonly examined 
but this has been observed in an unrelated material [l2T | 
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and without comment in a similar material ||- These 
features will be explained here. 

Scthna et al. have argued that one should be able to 
trade realistic microscopic degrees of freedom for a set 
of rules as part of a renormalization procedure to study 
different classes of hysteretic behavior at long length and 
time scales [l3| . While some multiparameter micromag- 
netic (time-dependent Landau-Lifshitz-Gilbert equation) 
models can produce hysteresis loops similar to those ob- 
served here |l4j . we here present a simpler model which 
explains these features in a simple, physically transparent 
manner and, we believe, captures much of the essential 
physics of these samples. 

Stochastic cellular automata models, like we use here, 
have been used successfully before to model nucleation 
barriers 0. Ours involves spins - representing individ- 
ual magnetic bits that are either up or down - on a 2D 
triangular lattice, as these materials are designed to work 
for 2D storage. The key model ingredients are as follows: 

1. Anisotropy field: Each spin is initially assigned an 
anisotropy field, defined as Hf. = 2K/((j,qM s ) where K 
is the anisotropy energy. The field — Hk can be thought 
of as the field at which an isolated spin would flip from 
up to down, and conversely for -\-Hu- Thus each spin 
has a coercivity and is inherently hysteretic. We assign 
different Hk values to each spin, typically according to a 
lognormal distribution. 

2. Energy barrier: At each step, we first treat each 
spin independently and compute an energy barrier for it 
to flip down. For an up-spin, in the case of magnetic 
field H > -H k , we use AE = ^(1 + -^-) 2 , where A is 
a constant, and for a down-spin in the case of H < Hk 
there is simply a sign change inside the quadratic [l5j| . 
The barrier is otherwise. 

3. Attempt probability: We evaluate the probabil- 
ity of the individual spin attempting to flip, P = 
exp(—AE/kBT), where T is an input parameter, and 
determine whether the spin attempts to flip. We do this 
for every spin on the lattice and keep a list of those that 
attempt. 

4. Cooperative rules: We then apply rules which re- 
quire cooperativity between spins and establish a nucle- 
ation barrier. The model we use is as follows: For a spin 
to flip from up to down, it must have at least three neigh- 
bors that are already down or will flip down in the same 
step (and similarly going from down to up). Thus the 
smallest size group of spins that can flip down is 7, as in 
Fig. [21 This effectively builds in a nucleation barrier and 
allows us to explore what effect this has on the hystere- 
sis loop. These rules play a diminishing role after initial 
nucleation. 

5. Magnetization and field: After updating every site 
on the lattice, we calculate and output the magnetization 
and step the magnetic field. This represents one simula- 
tion step. We then repeat this procedure until reaching 
saturation. 
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FIG. 2: Left, a group of seven spins as shown, where each of 
the seven attempts (small arrows) to flip down in the same 
step, is the smallest group of spins that can stably flip down 
(right) in our rules. Thus our rules provide a nucleation bar- 
rier to formation of a domain of down spins near M=l, and 
similarly when near M=-l for formation of a domain of up 
spins. 



It is important to note that we do not wait until 
the sample equilibrates before stepping the field. We 
thus link a simulation step with some small amount of 
real time, t, which has been done successfully before to 
model kinetic effects in hysteresis loops 16] . The rele- 
vant timescale for this may be the "attempt frequency" 
for crossing an energy barrier, which is generally consid- 
ered to be on the order of a nanosecond |l7j . 

At the most basic level, we show that without coop- 
erativity, our rules yield hysteresis loops that are quite 
broad. These loops acquire a sudden drop in magneti- 
zation due to the nucleation barrier when we apply the 
cooperative rules (Fig. Efa)). We also find that the tem- 
perature dependence of this effect is somewhat unusual 
(Fig-EHb)), in that the drop is actually sharper at higher 
temperature and both the drop and the tail are broader 
at low temperature. 

The explanation for these effects follows: The sudden 
drop in the hysteresis loop is due to the relief of the con- 
straints of the cooperative rules following nucleation of 
domains. That is, once the nucleation barrier is crossed, 
domains can grow, and it is much easier to add additional 
spins to a domain (a spin that is still up neighboring a 
down domain already has some down neighbors and so 
it is easier for it to have sufficient down neighbors to flip 
down) than it was to nucleate the domain in the first 
place. This is similar to a previous description of re- 
versal in terms of two coercivities, a nucleation field H n 
necessary to nucleate reversed domains and a field H p 
necessary to overcome domain wall pinning. If H n > H p 
reversal can occur very rapidly following nucleation 0,0] • 
However, our model is different in that we do not have 
a single pinning coercivity. We refer to the large cas- 
cade of spin-flips following nucleation as "avalanching," 
as described in random-field Ising models |l8j . although 
there is a distinction in that in the Ising model, reversed 
domains tend to pull their neighbors along with them, 
while here the only driving force is the external field. 

This avalanching applies to the temperature depen- 



3 




FIG. 3: (a) Descending half of hysteresis loops with and with- 
out cooperative rules for two different temperatures, T=0.04 
(thicker lines) and T=0.1 (thinner lines). The portion for in- 
creasing H is symmetric. The broad curves (solid lines) are the 
results when all the spins act independently; the sharp curves 
(dotted lines) show the effect of the nucleation barrier estab- 
lished by our cooperative rules, (b) A sequence of descending 
half hysteresis loops at different temperatures (T=0.03, 0.1, 
and 0.3 from left to right). Notice the overall sharpening at 
higher T and the much broader tail at lower T. (c) Calculated 
first-order reversal curves at T=0.1. Shown here are reversals 
at M=0.9, 0.5 and 0, but the curves essentially all overlap 
because in all three, even following reversal, magnetization 
continues to drop nearly to M=-l before stabilizing, (d) For 
T=0.05, first-order reversal curves are more reasonable. Note 
that the magnetization does drop following reversal. We find 
that the amount of this drop decreases at lower temperatures 
or for more broad distributions of Hk ■ 

dence of the drop in that at high temperature, more 
thermal energy is available to allow the avalanches to 
proceed. More thermal energy means bigger avalanches 
and sharper hysteresis loops, because it helps to over- 
come outlying Hk values that would tend to pin do- 
main walls. This pinning is possible because the rules 
do provide some constraint even after nucleation. Higher 
temperature means the disorder is less important and 
so avalanches are bigger. Low internal disorder has 
been shown to produce sharp hysteresis loops with large 
avalanches before 0, 0] , but our model includes both 
internal disorder and the effects of thermal fluctuations, 
and shows how temperature affects this sharpness. 

The broadening of the tail at low temperatures arises 
from a freezing in of the hysteresis loop to follow the ac- 
tual distribution of Hk at low temperatures. That is, in 
the tail of the hysteresis loop, our rules play little role 
(most remaining up spins have enough down neighbors 
they can flip as soon as they "want") and spins act in- 
dependently, so the tail samples the distribution of Hk 
and it must follow it more closely at low temperatures. 
This can be seen in Fig.02by the rejoining of the T=0.04 
curves in the tail region. This tail behavior emerges natu- 



rally from our model, ft has been observed before and its 
origin has been hinted at in phenomenological Preisach 
modeling || but we have provided a simple physical ex- 
planation. 

We also find (Fig. 03c), (d)) that our model does pro- 
duce a drop in the first-order reversal curves. At rea- 
sonably high temperatures (i.e. T=0.1 in our model) it 
greatly exaggerates this effect: reversing the magnetic 
field at M=0.9, the magnetization continues to fall al- 
most to M=-l before flattening out (Fig. Efc)). Lower 
temperatures or broader distributions of Hk make this ef- 
fect more reasonable (Fig. 02c)). Additionally, our model 
predicts that we can make this effect smaller by increas- 
ing the field sweep rate. That is, taking bigger steps in 
magnetic field gets back to the field necessary to stop the 
reversal curves from dropping sooner, so they do not fall 
as far. It is simply due to avalanching once a domain 
nucleates, it can grow easily and it requires getting back 
to a significantly less negative field to stop the avalanche 
of domain growth. Consider, for example, a spin with a 
low Hk which, if isolated, would be easy to flip. If it has 
several higher Hk neighbors it will not flip until long after 
it would like to flip. But once its neighbors flip it is quite 
happy to flip, even if the magnetic field reversal begins, 
since it is long past the point at which it would have liked 
to flip on its own. This also means that in our model, this 
asymmetric tail is present even when the distribution of 
anisotropy fields is symmetric (i.e. Gaussian). 

With this explanation in mind, experimentally, one 
should be able to allow the sample to equilibrate before 
taking a reversal curve and the effect should disappear, 
since it is a kinetic effect. On the other hand, slowing 
the field sweep rate along the reversal curve should allow 
avalanches to proceed further and the effect should get 
bigger. This was tested experimentally and confirmed for 
sample (a) (Fig. 0} and the drop also appears for sam- 
ple (b) if the field sweep rate is slowed enough j^jj. If 
we assume our theoretical model is capturing this effect 
accurately, r provides the conversion between simulation 
steps and experimental time. Thus a crude estimate of 
the effective r for these samples is on the order of a few 
milliseconds. This is obviously much larger than the typi- 
cal value of around a nanosecond and would suggest that 
these samples are unusual. This is consistent with the 
observation that this drop is not typically observed in 
experimental reversal curves. 

To check our prediction of tail broadening at low tem- 
perature, we used a Quantum Design Superconducting 
Quantum Interference Device (SQUID) magnetometer in 
the standard (DC) operating configuration with the sam- 
ple perpendicular to the applied field. Our procedure 
was to measure a full hysteresis loop at 5K and then just 
the descending half of the loop at 300K, to compare the 
broadening of the tail. Results are shown in Fig 03 both 
samples show a broadening of the tail at lower temper- 
atures in qualitative agreement with the predictions of 
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FIG. 4: Close-ups of the region where the reversal curves 
drop show the time-dependence of the effect for the sample 
of Fig^a). Essentially, when there is more pausing before 
or during measurement, the magnetization drops farther, (a) 
The normal reversal curves are taken by decreasing the field 
in 0.1 second steps, then increasing and averaging for 0.13s at 
each field point, (b) Averaging for 3s at each field point gives 
avalanches more time to proceed and the magnetization drops 
further, (c) Equilibrating for 60s prior to taking data on the 
reversal curve eliminates the effect but puts the magnetization 
lower at a given field than (a) and (b). (d) As (b) but the 
field sweep time prior to reversal is also slowed to 3s. 
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FIG. 5: Descending half of the major hysteresis loops for our 
Co-Pt multilayer samples as measured at two different tem- 
peratures in a SQUID magnetometer. Curves with most neg- 
ative coercivity are 5 K; the other are at 300K. The 5K curves 
also show the ascending half of the major loop, (a) The same 
sample as Fig. [Ra,) shows a significantly larger coercivity at 
lower temperature (leftmost curve) and a broader tail, while 
sample (b) shows a slightly broader tail (leftmost curve) but 
little change in coercivity. 



Fig. EJb) . This is more pronounced for sample (a) . How- 
ever, our model's hysteresis loops have a coercivity that 
is strongly temperature-dependent, but this is true for 
only one of these two samples. 

Overall, our results are especially significant in that we 
qualitatively explain the experimentally observed drop in 
the first-order reversal curves as a kinetic effect due to 
avalanching and provide experimentally verified sugges- 
tions for reducing or increasing this effect. This drop 
has been observed in another material and attributed to 
thermal relaxation |12|- thermal growth of domains at 



fields below that necessary to induce domain depinning. 
While thermal energy allows this drop to be larger in 
our model, we would argue both on theoretical and ex- 
perimental grounds that our mechanism is distinct. Ex- 
perimentally, we find that kinetics strongly influence this 
effect. And in our model, even at T=0, there is still 
a nucleation barrier which can result in avalanches and 
produce a drop in reversal curves - the drop is simply 
smaller. 

Additionally, our cooperative rules provide hysteresis 
loops with an asymmetric tail. This asymmetry is due to 
nucleation, because the nucleation barrier has a profound 
effect on the initial drop in magnetization. The effects 
of our cooperative rules disappear by the tail, allowing 
the tail to still be broad and individual spins to act in- 
dependently. This effect has only been captured before 
in vastly more complicated models [l4| or by parameter- 
tuning |5j- 

There are two possible origins of the experimentally 
observed broad tail. The first is the mechanism we de- 
scribe here, which does not depend on having an asym- 
metric distribution of Hk- To our knowledge, this effect 
is absent in models dealing in terms of exchange cou- 
pling, as exchange coupling strong enough to produce 
a nucleation barrier will also tend to prevent a broad 
tail. However, strongly asymmetric distributions of Hk 
can produce an asymmetry without hysteresis coopera- 
tivity. Independent measurements of the distribution of 
anisotropy fields will be very useful to distinguish be- 
tween these two cases. 

Our model suggests the importance of inherently hys- 
teretic spins, or hysterons, rather than Ising spins, where 
hysteresis comes solely from interactions between spins. 
In samples like these, with strong uniaxial anisotropy, 
there may be two separate energy scales, one associated 
with the energy barrier to flip individual grains, and 
one governing the exchange coupling between neighbor- 
ing grains which determines a cooperative nucleation bar- 
rier. The fact that nucleation is an irreversible process, 
as described above, helps to create a broad tail in our 
model despite the absence of any significant rounding at 
the top, and without invoking asymmetric distributions. 
Preliminary work suggests that this may remain true in 
an interacting hysteron model. On the other hand, mod- 
els based on Ising spins may not capture this behavior, 
as the deviation from fully polarized state at the top and 
bottom of the hysteresis loop would be due to reversible 
flips of outliers and thus would be symmetric. Hence, we 
suggest that interacting hysteron models, or cooperative 
nucleation models of the type discussed here, will be bet- 
ter for describing the hysteresis in materials with broad 
asymmetric tails. 

In conclusion, our model allows us to qualitatively ex- 
plain important features in the hysteresis loops of these 
samples and predict kinetic effects observed experimen- 
tally. We explain the drop in the reversal curves and the 
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broad tail. Our results suggest that exploring Ising-type 
models with exchange coupling and inherently hysteretic 
spins may give new insight into materials like these. 
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